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Abstract: We present the full details of a calculation at next-to- leading order of the momen- 
tum diffusion coefficient of a heavy quark in a hot, weakly coupled, QCD plasma. Corrections 
arise at 0{gs); physically they represent interference between overlapping scatterings, as well 
as soft, electric scale {p ~ gT) gauge field physics, which we treat using the hard thermal loop 
(HTL) effective theory. In 3-color, 3-flavor QCD, the momentum diffusion constant of a funda- 
mental representation heavy quark at NLO is k = ^a^T^{ln 0.07 428 +1.9026 gs). We ex- 
tend the computation to a heavy fundamental representation "probe" quark in large Nc , AA=4 
Super Yang-Mills theory, where the result is = (\n + 0.4304 + 0.8010\/A^ 

(where A = i^^A^c is the t'Hooft coupling). In the absence of some resummation technique, 
the convergence of perturbation theory is poor. 
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1. Introduction 



In the earliest stages of the Big Bang the universe was a relativistic plasma. We can now 
produce such a relativistic plasma in the lab via heavy ion collisions. In both cases the plasma 
is transient -in the early universe it lasted only around 10^® seconds and in a heavy ion collision 
it lasts little more than 10~^^ seconds. Since a system which remains always in equilibrium 
leaves essentially no traces of its earlier state, the most interesting physics in both situations 
is nonequilibrium physics. And in the early universe at temperatures above lO's of GeV 
[relevant for electroweak baryogenesis [1], leptogenesis [2], gravitino production [3], moduli 
production and destruction, and other relics] that plasma was weakly coupled. The same 
is true in principle for the early development of extremely high energy heavy ion collisions, 
though it is an open question whether this is a reasonable treatment at available energies. 

The natural language to study such plasmas is nonequilibrium quantum field theory. For 
plasmas near equilibrium (relevant for most of these interesting problems) one can study lin- 
ear deviations from (local) equilibrium by studying unequal-time equilibrium correlations and 
using linear response theory [4]. The tools for defining equilibrium finite temperature field 
theory and for performing weak coupling expansions have been known for over 40 years [5,6]. 
Nevertheless, our ability to calculate thermal and nonequilibrium phenomena is surprisingly 
immature. For thermodynamical properties we now understand how to compute pertur- 
batively very well; for instance, the best known quantity, the thermodynamic pressure, is 
known past fifth order [7-9]. We also understand that to compute real-time processes we 
need to perform a resummation of certain plasma effects via the so-called Hard Thermal 
Loops (HTL's) [10]. 

However surprisingly few gauge invariant real-time correlation functions have been com- 
puted even at leading order, and even fewer are known beyond this level. In particular, 
transport coefficients - shear viscosity, baryon number diffusion, electrical conductivity, heavy 
quark diffusion, bulk viscosity, and so on - are of considerable importance, since they describe 
the relaxation of a system which is relatively close to equilibrium. However, even leading or- 
der calculations of these quantities only became available quite recently, and none of them 
are known beyond leading order. 

This is a major gap in our understanding of finite temperature field theory. In particular, 
it has been known for some time that the rate of convergence of the perturbative series in 
the thermal theory can be much worse than in the vacuum theory. For instance, Braaten and 
Nieto [8] argued that the convergence of perturbation theory was poor unless < 0.1, due 
to physics at the so-called electric screening scale ~ gT (where the loop expansion really is 
only an expansion into powers of g.) It may be possible to rescue, or at least improve, the 
convergence of perturbation theory using various resummation techniques [11, 12]. However 
neither of these issues has been explored for dynamical (unequal time) quantities. In this 
context exactly one interesting gauge invariant quantity is known at next-to-leading order; 
the deeply virtual dilepton production rate [13]. However this quantity involves short physical 
time scales and is not sensitive to the electric gT scale, and is therefore not representative of 



- 1 - 



the other transport coefficients, which are. For these "soft-sensitive" quantities, we do not 
know what a next-to-leading order computation involves, we do not know how well the series 
will converge, and we do not know what resummation techniques might be available or how 
well they may work. 

The simplest such quantity is the heavy quark diffusion coefficient, first computed at 
leading order in [14]. This is a quantity of phenomenological interest, which sets the rate at 
which the velocity of a heavy quark equilibrates with that of its environment. In a recent Let- 
ter [15] we have presented the result of a computation of the heavy quark diffusion coefficient 
in full QCD to next-to-leading order. Here we present the full details of this calculation, as 
well as its extension to AA=4 super- Yang-Mills theory (AA=4 SYM). An outline of the paper 
is as follows. In the next section we summarize the problem, discuss the relevant physics, 
and present the results. Then Section |^ presents our approach (HTL effective field theory) 
and reviews the leading order calculation. Section ^ presents the body of the calculation of 
the NLO diffusion coefficient, and Section |5| extends this treatment to AA=4 Super- Yang- Mills 
theory. A technical appendix discusses the analytical properties of some of the functions we 
encounter in the calculation. 

2. Overview and results 

2.1 Definition of heavy quark diffusion 

A heavy quark, M ^ T, in or near equilibrium has a typical momentum squared ~ MT ^ 
large compared to the plasma scale and it therefore takes a parametrically long time for 
the momentum to change appreciably. This means that momentum changes accumulate from 
many uncorrelated "kicks," so on long time scales p will evolve via Langevin dynamics, 

^ = -riDPi + m , imC.it')) = K5,,5{t - t') . (2.1) 

The relaxation rate r]£) and the momentum diffusion coefficient k are related by a fluctuation- 
dissipation relation, r]D = jMT^ which follows on general thermodynamical grounds. Thus 
the dynamics of the nonrelativistic heavy quark is completely set by the single parameter k, 
which we compute to next-to-leading order. 

2.2 Qualitative origin of NLO effects 

At leading order in the weak coupling expansion, the momentum diffusion coefficient is set by 
the t-channel Coulomb scattering processes illustrated in Fig. 0, in which the scattering target 
can be a light quark or a gluon (Compton-like processes are suppressed for nonrelativistic 
heavy quarks.) An important feature of weakly coupled plasmas, relativistic and nonrela- 
tivistic, is that the total rate for Coulomb scattering is quadratically divergent in the limit of 
small momentum transfer |q| [16]. Such a divergence is of course unphysical in a medium of 
charged particles, and is cut off by screening effects^ at the momentum scale q ~ rrin ~ gT. 

^In relativistic plasmas there still remains a logarithmic divergence due to un-screened magnetic scatterings, 
but this cancels in physical quantities and will not concern us here. 
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Figure 1: Processes responsible for heavy quark diffusion at leading order: Coulombic scattering 
between a heavy quark (double line) and either a quark or a gluon. 



The momentum diffusion coefficient k is sensitive not to the total scattering rate, but to 
the weighted average, 

2.2 



- / d3,.^r(g) , 



where dT{q)/d^q denotes the differential probability per unit time for the momentum of the 
heavy quark to change by q. The two additional powers of q"^ present in Eq. (^]^) reduce 
the quadratic divergence of the total rate / d?qdV/d?q to a logarithmic divergence, cut off 
at g ~ mj3. The logarithm reflects the fact that the heavy quark momentum diffuses due 
to a range of momentum transfers, from the many soft scattering events, which individually 
have q ~ gT but occur on a rate Fsoft ~ g^T, to the rare hard scattering events with q ^ T, 
occurring on a rate Fhard ~ d'^T. 

Usually interaction corrections involve powers of g'^, and for large momentum transfer 
processes this is true. But the presence of the Coulombic divergence, cut off by screening 
effects, means that the details of soft momentum exchange and plasma screening are relevant 
to heavy quark diffusion already at leading order. Now consider the contribution from gluons 
to the Debye screening scale [17]: 



d^k nB{k) 
(27r)3 k 



+ (fermion contribution) . 



(2.3) 



Here = A'^c = 3 is the trace normalization of the adjoint representation and is the Bose 
distribution function. At small k the integral behaves like ~ g^T J d\k\, due to the singular 
nature of the Bose-Einstein function nB{k) ~ T/k. Thus the contribution from gluons with 
k ~ gT represents 0{g) of the total strength of plasma screening. 

However, the soft gluon contribution is not computed correctly by the above expression. 
The derivation of Eq. (2.3) assumed free massless propagation of the particles responsible 
for screening, but gluons with 0{gT) momenta themselves experience 0(1) screening effects. 
Therefore one must recompute that part of plasma screening which arises from the 0{gT) 
gluons. This calculation is complicated by the fact that interactions between soft gluons are 
also strongly modified by plasma effects, described by the HTL effective theory [10,18]. Thus, 
since a relative 0{g) fraction of Eq. (|2.3| ) arises from soft gluons, a correct treatment of the 
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effect they produce will produce an 0{g) correction to the physics which cuts off the infrared 
logarithm in /t, and therefore an 0{g) correction to the result. 

For similar reasons, in an 0{g) fraction of soft scattering events, the plasma particle 
which strikes the heavy quark is itself a soft gluon with momentum ~ gT. But the dispersion 
and spectral weight of such gluons are strongly modified by the plasma and this contribution 
to the "target density" must also be reconsidered. 



K p 




Diagram 4 Diagram 5 Diagram 6 Diagram 7 



Figure 2: Diagrams which can interfere when two scattering processes overlap. Here K and P are 
the initial momenta of two distinct hard particles, which can be quarks or gluons. 

Other NLO corrections can be expected to arise from interference between overlapping 
scattering events, as illustrated in Fig. |2|. Since the total scattering rate is ~ g'^T, both for 
the heavy quark and for the light particles which scatter it, and since small angle scatterings 
can have a duration up to ~ ^/gT (this can be read off, for instance, from the virtuality of the 
exchanged gluon), an 0{g) fraction of scattering events overlap with another scattering. Thus 
the individual diagrams in Fig. |2| are only down by 0{g) relative to the leading processes. 
Corresponding virtual corrections to the leading order scattering process will naturally also 
arise at this order. However, the appropriate question to ask is not about the probability 
for scattering events to overlap, but rather how much do they interfere with each other. In 
QED, the small angle scattering of a particle occurring during another scattering event has 
little impact on that event, and correspondingly one finds a parametric cancellation between 
diagrams 1 and 2 of Fig. |2| (and between diagrams 4 and 5 and between diagrams 6 and 7, and 
among the associated virtual processes; diagram 3 is absent altogether in QED.) However, 
in QCD, instead of a cancellation one gets a commutator of group theory factors: the very 
frequent soft scatterings which occur in the plasma do matter, because they change the colors 
of the particles. 

Note however that this source of 0{g) NLO corrections and the preceding one are not 
clearly distinct. Indeed, diagram 3 of Fig. § can be understood as the special case of the 
diagram in Fig. |l], where the external gluons are soft and in the Landau cut. This suggests 
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that, to make the quahtative discussion here more precise, we will need to perform a careful 
diagrammatic approach based on power counting. There is one common feature of the sources 
for correction we have listed, however; all involve the influence of soft gluons. This observation 
suggests that the calculation may be rephrased in terms of an effective theory of gT scale 
physics, in which the hard scale ~ T has been integrated out. This is precisely Braaten and 
Pisarski's HTL effective theory [10]. Carrying out a careful diagrammatic calculation within 
this effective theory is the subject of the body of this paper; in the remainder of this section 
we will present the results. 

2.3 Results: QCD 

The squared matrix elements for the processes of Fig. ||, summed over the initial and final 
states of the light scattering targets and final states of the heavy quark, and averaged over 
the initial states of the heavy quark, have been evaluated in [19], yielding 



+NfnF{k){l-nF{k)) (2-^) , 



(2.4) 



where C^^ = | in QCD is the quadratic Casimir of the heavy quark representation, and 
rrio = \/1.5gT in QCD with A'^f=3 flavors of light quarks. Formally taking m^, <C T, the 
integral is dominated by /c ~ T and q in the logarithmic range nio ^ q ^ T. The leading 
behavior in g of Eq. ( |2.4| ) can be obtained from the leading behavior in mj^/k'^ of the q 
integral. Making room for the next-to-leading order correction C, the result can be written: 

C,g'T^ 



K- 



ISvr 



Nc+— 
2 



, 2r ^ 

In +^ 



A^fln2 Nf.mn ^ oxi , ^ 

+^^+^Y^C + 0{g')] . (2.5) 



Here ^ = ^ — 7^ -|- 1^ ^ —0.64718. The leading order part of Eq. ( |2.5| ) was given explicitly in 
[19] (it could also have been extracted from the nonrelativistic limit of earlier results [14,20].) 
The dependence of the next-to-leading order correction on physical parameters is contained in 
the coefficient multiplying C, which itself is a pure number: all of the above-mentioned next- 
to-leading order corrections depend on physical parameters in the same way as an 0{m]j/T) 
fraction of the gluon contribution to k^*^. 

Expression Eq. ( |2.4D itself contains 0{g) corrections, giving rise to a rather trivial con- 
tribution^ to C, 62^2 = ^ 0.8356. It arises wholly from the k ~ gT region of the gluon 
contribution to Eq. ( |2.4|) , where the result of the q integration is poorly described by the 
leading term of its m^/k'^ expansion, which was used to obtain the leading order behavior 
Eq. ( |2.5D . Although slightly tedious, the evaluation of is entirely straightforward and we 
do not present it here. In section ^ we compute the difference between the full next-to-leading 
order momentum diffusion coefficient, and what is already incorporated in , and obtain 
C ~ 1.4946. Thus C = C2-.2 + C ^ 2.3302. 



^In [15] this contribution was named Cg^ . 
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Figure 3: Comparison of leading and NLO results for iVf — 3 QCD as a function of coupling. 

Our result Eq. ( p.5| ) is plotted in Fig. |3[ A simple-minded estimate of the regime of 
validity of perturbation theory can be given by equating the size of the correction to the size 
of the leading-order result. What is usually referred to in the literature as being the leading 
order result is Eq. ( |2.4D , numerically integrated at a given value of the coupling (this is the 
curve called "leading order" in Fig. ^): the correction becomes as large as this leading order 
result when ^ 0.04. This suggests that at that point perturbation theory starts to get into 
trouble. For this reason, and as should be clearly suggested by the plot, we do not believe 
that our calculation can be directly used as an "improvement" to the determination of k in 
the context of heavy ion collisions, where phenomenologically realistic values of the coupling 
are in the range as ~ 0.3 — 0.5. Rather our results signal difficulties with the approach. 

Nevertheless we would not like to sound overly pessimistic and conclude that our results 
signal that no prediction beyond as = 0.05 is possible. Rather, the real question now is how 
large the higher order corrections are, and more pertinently, which parts of C may duplicate 
themselves in higher-order terms, in some more or less predictable (and therefore resummable) 
fashion. 

Consider for instance the difference between the two lowest curves of Fig. |3|, which is 
attributable to C2—,2, up to terms which are of yet higher order in the rrio/T expansion of 
Eq. ( |2.4| ). This contribution, which can be evaluated knowing only the tree-level matrix 
elements with massless external states (and HTL corrections resummed on the exchanged 
gluon), is better described as an "ambiguity" in the leading-order result rather than as a 
correction to it. This ambiguity is large because the Coulomb scattering processes against 
soft gluons (which give the small k contribution to Eq. ( |2.4| )) are poorly described by the 
leading term of an rrin/T expansion. This is unrelated to the question of whether these 
processes are correctly described by the right-hand side of which is the most pertinent 
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question to ask if we are concerned with higher order terms. Actuahy, our calculation tells 
us that the effect of the HTL changes in the dispersion relation and interaction strength of 
the soft gluons is relatively modest, essentially given by the pole-pole contribution of section 



4.3.4 of order C(^)^poie-poie ~ —0.20. Thus it appears that this region of phase space is not so 
poorly described by Eq. (2.4), and that simply defining this full expression to be the leading 
order result should provide a reasonable resummation of the contribution C2^2- 

Along the same lines, there is another contribution to C which has a simple physical 
interpretation and which would be easy to include into the leading order calculation. As 



we will discuss in Subsection |4.2| , just over half of the remaining correction arises from a 
shift in the (real) Debye screened propagator l/(g^ + m^) appearing in Eq. (|2.4| ). This can 
be understood as an NLO momentum dependence in the Debye screening mass m^. This 
momentum dependence can be resummed in a few ways. One way would be to solve for it 
in the 3D Euclidean effective theory nonperturbatively or via a gap equation (though this 
approach appears to be special to heavy quark diffusion, where the exchange momentum is 
strictly spatial). Another method would be to replace in Eq. (2.4) with the full momentum- 



dependent leading-order self-energy at general q (though this procedure does not appear to 
be gauge invariant). 

For these reasons we think with some optimism that two thirds of the difference between 
the lowest and highest curves in Fig. ^ can be absorbed in a relatively simple systematic 
resummation scheme, and at most only one third represents complicated physics that will be 
really difficult to resum. Such a resummation scheme might then extend the reach of pertur- 
bation theory to Os ~ 0.15 or so-high enough for almost all cosmological applications, though 
still not enough to be much use for heavy ion physics. Clearly such issues of resummation 
are an interesting problem for future work. 

2.4 AA=4 super Yang-Mills 

Results for heavy quark diffusion in the strong coupling regime of J\f=4: super Yang-Mills 
(SYM) have been obtained in the literature, exploiting the AdS/CFT correspondence [21,22]. 
It seems interesting to study the heavy quark momentum diffusion coefficient also at weak 
coupling in this theory, for the purpose of comparison between the two regimes and for 
comparison between this theory and ordinary QCD. The leading order result at weak coupling 
has previously been given [23], and here we present the next-to-leading order correction it 
receives.'^ 

We begin with a brief description of AA=4 SYM and of heavy quarks in this theory. In 
addition to the gauge field A^, N=^ SYM contains four Weyl fermions and six real scalars, all 
transforming under the adjoint representation of the gauge group. The theory contains a single 
dimensionless coupling constant (7, which sets the strength of the gauge, Yukawa and scalar 
interactions; the Lagrangian is completely determined by the supersymmetry [24]. The strong 



^We present our results at leading order in the large A^c expansion. Strictly speaking, the theory at finite 
= \/Nc and with added fundamental matter has a Landau pole; however it is valid perturbatively. To 



generalize Eq. (2.6) and Eq. (2.7) to finite A'^c, multiply the righthand sides by 2Ch / Nc 
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Figure 4: Leading order diagrams involving scalars, which are present in Af—4: SYM but not in 



coupling results are obtained in the large A'c limit of the theory with gauge group SU(A'^c)i 
and for this reason we will express our results in terms of the t'Hooft coupling A = g^Nc- 
What is meant by a "heavy quark" in this theory is an J\f=2 massive hypermultiplet added 
to it, transforming under the fundamental representation of the gauge group. In terms of 
the M = 2 field content of M=^ SYM, this heavy hypermultiplet is minimally coupled to 
the M = 2 gauge multiplet of the theory, but is not directly coupled to the massless matter 
hypermultiplet. This is the conventional setup employed in AdS/CFT studies [21]. 

In the large mass M ^ T limit, processes which would change the identity of the heavy 
particles, from heavy quarks to heavy scalars and vice- versa, are suppressed [23]. Thus the 
heavy fermion carries an approximately conserved U(l) charge and it makes sense to speak 
about its momentum diffusion coefficient, with no reference to its scalar superpartners. In 
addition to the scattering processes depicted in Fig. in this theory at leading order there 
are scattering processes involving light scalars, depicted in Fig. ^ Including these processes, 
the leading order momentum diffusion coefficient k can be written [23]: 



where = 2XT and m| = AT. The first line describes Coulomb scattering against gluons, 
the second line describes Coulomb scattering against five real scalars and Coulomb plus 



QCD. 




2 



(2.6) 



8 



0.15 



0.1 



CM 
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Next-to-leading order - 
Leading order - 
Truncated leading order 



Figure 5: Comparison of next-to-leading and leading order results for heavy quark diffusion in 7V=4 
SYM theory. 



Yukawa-Compton scatterings against one real scalar^. The third line is Coulomb scattering 
of fermions, the fourth line is the scalar-mediated conversion of light gluons to light scalars and 
vice-versa, and the last line contains the scalar-mediated scatterings against light fermions. 
The integrals of Eq. (|2.6| ) were evaluated to leading order in rrio/T ~ \/A in [23]; making 
room for the next-to-leading order contribution the result can be written: 



X^T^ 2T \ 1 , ^ 

In h h - In 2 

Gvr \ rrio 2 3 



(2.7) 



with ^ as defined below Eq. (2.5). As for QCD, a rather trivial contribution 



(sym) 



15 
2-K 



— ^ ~ 1.7121 to C^^™) arises from the expansion of Eq. ( |2.6| ) to next-to-leading order in 

coming from the k ~ ^/XT region of processes involving external bosons. Another part 

of C is precisely the same as the "difficult" part of the QCD calculation, C'('3c«) ~ 1.4946. 

The remainder C^^^^^) ~ 0.19172 is calculated in section |. Thus C^^^*^) = c':fZ2^ + C^'^^^) + 
C(sYM) ^ 3,3984. 

This result is plotted in Fig. ^. As the figure shows, the convergence of the series is 
somewhat better than in QCD. Viewed as a function of the t'Hooft coupling A the correction 
is 100% of the leading value for A ~ 1.71, to be compared with the QCD value of as = 0.033, 
which is A = 1.26. Viewed as a function of the Debye screening scale, the comparative 
convergence in SYM would seem even better; the 100% correction point occurs for mj, = 



''Apparently the Yukawa-Compton scattering processes against scalars (the second diagram in Fig. ^ were 
not included in the calculation of Vuorinen and Chesler [23]. This caused an error in their determination of 



the constant term in Eq. (2.7); they found 



rather than ^. 
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1.85T, whereas it is mo = y I^T = 0.79T. However, estimates of this kind can sometimes 
be deceptive: the fact that at fixed rrio/T the NLO correction in SYM is comparatively 
smaller than in QCD can be mostly attributed to the leading order result being stronger 
in SYM, due to its larger number of matter fields (which also scatter with a larger group 
theory factor, being in the adjoint representation.) However, this additional physics that is 
present in SYM suffers from relatively modest NLO corrections, the most severe corrections 
still being associated with soft gluons. Thus, although the leading order result has a wider 
range of validity in SYM than in QCD, it seems reasonable to expect the range of validity of 
the NLO correction itself to be roughly similar in SYM and QCD, in terms of m^/T. 



3. Effective theory and leading-order analysis 

Our first step towards a rigorous analysis of the momentum diffusion coefficient k is a non- 
perturbative definition [22] in terms of a force-force (electric field-electric field) correlator 
with Wilson lines connecting the electric fields: 

K = ^ dt l^H{W{t; -oo)t Et{t)t% W{t- 0) (0)4 W{Q- -cx))) . (3.1) 

The trace runs over the representation of the gauge group of the heavy quark, and the Wilson 
lines W act on this representation. Intuitively, Eq. (|3.1| ) is exactly the force-force correlator 
of Eq. ( |2.lD , with the forces given by electric fields and the Wilson line representing the gauge 
rotation of the heavy quark due to propagation, which ensures gauge invariance. Because of 
operator ordering issues, the Wilson lines shown are not equivalent to connecting the E fields 
with an adjoint Wilson line. The Wilson lines also incorporate the effect of the heavy charge 
on the plasma (which is why they must go back to time — oo). In general they introduce 
nontrivial representation dependence into the heavy quark diffusion constant, and in fact 
such Wilson lines are even required in QED (diffusion of ions in a QED plasma depends on 
the ionic charge Z in a more complicated way than only because of these Wilson lines, 
which account for the reaction of the plasma to the presence of the charge). However we will 
see that to the order we work here, they can be replaced by an adjoint Wilson line. 

Our approach is to calculate this correlation function within HTL effective field theory. 
This is an infrared effective description valid below a cutoff scale ~ T which describes gauge 
fields and fermion fields, resumming into the Lagrangian certain 0{T'^) plasma effects. Per- 
turbatively this introduces corrections to the propagators and vertices which become 0(1) 
at the scale ~ gT, which is a natural scale in HTL effective theory. The effective theory 
requires matching to the full thermal theory in the UV, requiring counterterms both in the 
Lagrangian parameters and in correlation functions such as Eq. ( |3.1| ). Perturbation theory 
within the HTL effective theory is expected to converge in powers of 5, which intuitively can 
be though of as the usual factor times a Bose statistical factor evaluated at the scale gT, 
UsigT) 1/g. The HTL calculation can also encounter infrared divergences, arising from 
the unscreened low-frequency magnetic gluons. The appearance of such an IR divergence 
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signals the breakdown of perturbation theory and the need for nonperturbative information 
about the ultrasoft magnetic sector. We expect such IR divergences at some finite order in 
perturbation theory, but this proves to be beyond the NLO level we consider here.^ 

Let us proceed with the leading order calculation of the correlator in Eq. (^T^). At leading 
order one may replace the Wilson lines with identity operators and use the noninteracting 
form of the electric field correlator,^ Ei = OiAq — d^Ai. The time integral means we need the 
result at zero frequency and so the d^Ai piece does not contribute^. Fourier transforming to 
the momentum basis, we rather immediately obtain 



Ma J (27r) 



r(7'(^Oa(u; = 0,(?)Aoa(0)) 



r9'G>o(0,(?) 



(3.2) 



The Wightman propagator is to be evaluated within the HTL effective theory. Using the KMS 
condition we can express the Wightman propagator in terms of the retarded propagator. 



(3.3) 



which is given in the HTL effective theory, in strict Coulomb gauge (which we use throughout), 

by 



GliP) 



with 



n?°(p) 



ml, 



ml, 



^In 
2 



l + r]\ 



+ 



1 — rj 
r/(l -r/2) 



p2 + n00(P) 
—i 



In 



1 + r/l \ i7rr;(l — r/^) 



1 — r/l 



0(1 -v') 



where t] = /p. Therefore the leading order momentum diffusion coefficient is 



lim 



dh 



,2T 



g'^Cam^ 



(27r)3^ uj 2g(g2+m2)2 



67r 



q^dq 



(g2 + 777,2)2 



(3.4) 



(3.5) 



(3.6) 



This integral is UV log divergent, indicating the need to perform a matching calculation. 
This is done by finding the result in the full theory, which gives Eq. (|2.4|) (without the m^ 



^We believe that most transport coefficients are sensitive to nonperturbative magnetic physics at 0{g'^), 
which is the relative contribution of these magnetic fields to transverse momentum diffusion for a moving 
particle {v ~ 1) in the soft electromagnetic fields of the HTL effective theory. Because the heavy quark 
considered here has v ^ 1, we believe magnetic physics arises at a higher order than 0{g^). 

®We work in [ — h++] metric contentions and 4-vector potential A'^ = (A", A') with A" the usual scalar 
and the usual vector potentials. 

'^This does not apply in gauges such as temporal gauge where the zero frequency gauge boson propagator 
can display singularities. 
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terms in the denominator). In the range gT <^ q <^T the two calculations agree; performing 
the k integration in Eq. ( |2.4| ) treating q <^ k reproduces Eq. (|3.6| ). The matching should be 
performed in some way so that the UV region is equivalent to Eq. ( |2.4D and the IR region is 
equivalent to Eq. ( |3.6D . One way of doing this would be to compute both Eq. ( |2.4| ) (without 
factors in denominators) and Eq. (|3.6| ) each in dimensional regularization and add them; 
the 1/e factors will cancel and the finite parts will give a consistent leading order result [25]. 
Alternately one can simply insert factors in the denominators in Eq. ( |2.4p (as has already 
been done) so that one expression is appropriate in the IR and UV. The error thus introduced 
at large q is only NNLO {0{g^)) and will not interfere with our NLO calculation. 



4. Details of the calculation 



We now proceed to push the leading order calculation of the last section to the next order in 
HTL perturbation theory. 

4.1 Formalism and diagrams 

The real-time correlator Eq. (3.1) can be expressed in terms of correlators of fields ordered 



along the Schwinger-Keldysh contour [6,26]. We find it convenient to use the Keldysh or ra 
basis (where one works in terms of the contour averaged or r and contour differenced or a 
fields rather than the fields on the upper 1 and lower 2 contours). Throughout we will be 
using a graphical notation for the propagators of this formalism, summarized in Table ||. We 
draw retarded (ra) and advanced (ar) propagators with an arrow on them, which points in 
the direction of the flow of time (thus towards the r index.) We draw the rr propagators 
with a double cut in the middle of it; there exists no aa propagator in this formalism. If one 
thinks of an rr propagator as carrying two outgoing arrows, leaving away from the cut, then 
an r field at a vertex will have an incoming arrow on it and an a field will have an outgoing 
arrow, leaving the vertex. Thus the ra assignments of the fields at the vertices can be readily 
recovered from our notation. The nonzero (tree- level) vertices in the Keldysh basis carry odd 
numbers of a indices, thus have an odd number of arrows leaving them (but arbitrarily many 
r indices, or incoming arrows.) Interaction vertices with one a index are precisely the same 
as the usual ones given in textbooks on quantum field theory at zero temperature, and those 
having three a indices are smaller by a factor |. This notation is the same as in [27]. 

Diagrammatic rules that generate the HTL effective vertices have recently been worked 
out for real time field theory in this basis [27]; we use the results and notation of that work, 
which are summarized in Fig. ^. In the HTL limit the hard degrees of freedom behave like 
classical point-like particles, which we draw as solid lines. The vertices (a) and (c) and 
the eikonal propagator (d) of the figure together describe the generation and propagation of 
disturbances of the hard particles' distribution functions in a background gauge field, and 
the two-point vertex (b) describes how these disturbances source gauge fields. These effects 
depend on the four-velocity v'^ = (1, v), = 1, of the hard particle, which has to be averaged 
over for every connected solid line that appears in a diagram; a factor of rnj^/T must also be 
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Symbol 



Notation 



GraiP) = Gr{P) 



Gar{P)=GA{P) 



Grr{P) 



Expression for a free scalar field 



p2 



p2 _^ jgpO 



Table 1: Graphical notation for real-time propagators in the Keldysh basis, and their expression for 
a free scalar field. In all cases the momentum P flows from right to left. 



(a) a mnr r 



(b) a^TT 



(c) a- 

(d) a- 

(e) r- 



/i, a 
b 



< 



P 



p 



—i 
V ■ P- 

2t:5{v ■ P) 



Figure 6: Feynman rules for the HTL theory, with ra indices explicitly shown: (a)-(c) give inter- 
action vertices and (d)-(e) give effective propagators for classical particles. All two-point functions 
are proportional to the identity in color space, (5°'', not explicitly shown. A factor {mjj/T) plus an 
integration J over the four-velocity must be assigned to every disjoint solid line appearing in 
a diagram. 



included. The rr propagator (e) describes statistical fluctuations in the number density of 
the charges, and enters precisely once in the calculation of HTL amplitudes with two external 
a gluons. The only HTL three point vertex that exists has the ra assignment shown in (c); 
there exist no HTL amplitudes with more than two external a gluons. Diagrams containing 
self-energy insertions on gluon propagators must be discarded, since we are already using 
the HTL-resummed gluon propagators Eq. (0). The apphcation of these rules reproduces 
calculations in classical Yang-Mills plasmas with point-like (nonabelian) charges; more details 
can be found in [27]. 

In the Keldysh basis, correlators involving Keldysh a fields with soft momenta p ~ gT 
are systematically down by powers of g, relative to similar correlators with Keldysh r fields 
(see e.g. [27]), implying that at NLO all fields entering Eq. ( |3.1| ) can be taken to be Keldysh 
r fields. This means that operator ordering issues actually are subleading, and that at NLO 
the Wilson lines in Eq. ( |3.lD can be traded for a single adjoint Wilson line. 

It is convenient to write the electric field operators as -E* = —d^A^ — DfA^. Using the 
equation of motion for a Wilson line, DtW{t] 0) = 0, we can then express 

W{t; 0)t {DtA^^it) tl) W{t; 0) = | 0)t A^-{t) W{t; 0)) , (4.1) 
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(A) 



(B) 



(C) 



(D) 



Figure 7: The Feynman diagrams that contribute to k^'"^, with assignments of momenta, Lorentz 
and color indices. All propagators are soft and HTL-resummed, and all interaction vertices include 
HTL corrections; the arrows denote only the momentum flow. The Lorentz indices of the gluons which 
connect to the heavy quark (shown as the double line) are all "0" . 



which contributes a total derivative to Eq. (|3.1| ) and can be dropped. Therefore the i?* can 
be replaced with — SM", and the desired correlation function becomes 



3dA 



dtid'A^^it) 



9^0^(0)) . 



ab 



(4.2) 



Expanding the Wilson line gives a series of correlators of ^4^^ fields. The diagrams we 
need at NLO are shown in Fig. |7[^ These diagrams are to be evaluated within the HTL 
effective theory, meaning that all propagators are HTL resummed and all vertices include 
HTL vertices-except for the vertices on the Wilson line, shown as the double line in the 
diagram. Naively there could be two more diagrams, corresponding to (A) and (C) but with 
fermions rather than gauge bosons in the loops; but these are suppressed by at least one 
factor of Uf/rLi, ~ g relative to the indicated diagrams and can be neglected at NLO. 

4.2 Real part of self-energy diagrams (A) and (C) 

Diagrams (A) and (C) of Fig. correspond to NLO self-energy corrections to the soft zero- 
frequency longitudinal gluon propagator, and can be decomposed into real and imaginary 
parts. We begin with the real part of the self-energy. Since it is needed only at zero frequency, 
it can be most convenient evaluated within the imaginary time formalism, in which the 
frequency integrals are replaced by discrete sums over the (imaginary) Matsubara frequencies 
UJn = 2TTnT, n integer [28]. In this formalism, no analytic continuation of any kind is required 
at zero external frequency, and we can directly analyze the discrete sum over the Matsubara 
frequencies. 

Because the HTL effective vertices vanish when all of their external frequencies are zero, 
and are subleading by (and at any rate, inappropriate) when one of their external momenta 



*This is true in any gauge in which the propagators show no divergent or pathological behavior in the 
zero-frequency limit. These includes the covariant or Coulomb gauges, but not the temporal axial gauge. 

^There is also a diagram which looks like (D) but with crossed gluon lines. It vanishes by rotational 
invariance. 
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carries a nonzero Matsubara frequency ^ T, the diagrams involving HTL vertices do not 
contribute at NLO. Similar cancellations occur for the transverse-transverse contribution with 
tree interaction vertices, because the relevant interaction vertex vanishes when all frequencies 
are zero, and the contribution of nonzero Matsubara frequencies only receives 0(g'^) correc- 
tions (the presence of the hard frequency scale in the loop propagators ensures that the 
self-energy corrections on the loop propagators are down by 0{g'^), and that the p depen- 
dence of the integral over spatial momenta can be expanded into integer powers of /T"^ 
when p <^T < uJn-) The diagrams with topology (C), involving four-point vertices, similarly 
do not contribute: the one with an HTL vertex is irrelevant at zero external frequency, and 
when the vertex is a tree vertex the propagator in the loop must be purely transverse (in the 
strict Coulomb gauge) since there is no interaction vertex involving only fields; but the 
HTL correction to this propagator at zero frequency vanishes. These simplifications, which 
are specific to the zero- frequency retarded self-energy (and to a lesser extent, to our use of 
strict Coulomb gauge), are perhaps best understood in terms of the dimensionally reduced 
effective theory [29]. 

We are thus left to evaluate the transverse-longitudinal loop with tree-level interaction 
vertices and HTL-resummed propagators. Only the contribution from the zero Matsubara 
frequency is needed. 



-g^N.T / 4p2 I 1 



2tt 



p^ 



m 



D 



1 



1 



{p ■ qf 

p2q2 J g2 ^2 _|_ j^2 



D 



P 



tan-1 ( ^ ) 



(4.3) 



where is shorthand for J ^^^^3 and where the arctangent takes values in [0, The same 
result could also have been obtained within the real-time formalism, albeit with somewhat 
more work. This 0{g) correction to the real part of the longitudinal gluon self-energy induces 
an 0{g) correction to k, 



Snne = ^J^p'U>'^'ip) 



1 



(p2 + m|, + (5n^(p))2 {p^ + m|, 



2 ^2 



(4.4) 



which, using the HTL approximation Il^^^{p) = m^vrT/p, yields a contribution to the di- 
mensionless coefficient C from Eq. (p.5|): 



C(A),Re = Svr 



P 



,(l+p2)3 



1 + 



P 



1 



■ sm 



-1 



P 



3 / tt \ 

-1 + -)=. 0.77198914., 



(4.5) 



where we have rescaled p to p/m^. Note that the integral is both IR and UV safe. 
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4.3 Self-energy (A): imaginary part 
4.3.1 Overview of the calculation 

The imaginary part of the gluon self-energy diagram (A) is probably the most technically 
challenging part of this calculation. Instead of calculating the {0{p^) term of the) imaginary 
part of the retarded self-energy, we find it more convenient to calculate directly the Wightman 
self-energy Il-^^^{P) = 2{l+nB{p^)) Imll^j^, which can be evaluated directly at zero frequency. 
There exists a finite temperature cutting rule, analogous to the familiar zero-temperature 
Cutkowski rule, which expresses this function in terms of a sum over diagrams that are 
divided into two parts by one cut [30,31]!°. rpj^g 

propagators traversed by the cut are to be 
evaluated as Wightman propagators, = {Gji{P) — G^(P))(1 ±n{p^)), for bosons and 

fermions respectively, and the "amplitudes" on each side of the cut are to be evaluated as the 
fully-retarded amplitudes of the real-time formalism, the retardation (e.g. time flow) being 
taken to be away from the cut, toward the external legs. In terms of the Keldysh ra basis, this 
means that all cut propagators attach to the neighboring vertices like Keldysh r fields, and 
external legs of the self-energy diagram should be considered as carrying Keldysh a indices. 
These fully-retarded amplitudes are the simplest analytic continuation of the imaginary-time 
amplitudes [32] (they are obtained by continuing all but one of the external momenta from the 
upper-half complex frequency plane.) A direct proof of the rule we use, within the real-time 
formalism, has also been given using the R/A formalism [33] (see also [34], section 3.6.); the 
cutting rule presented there is the same as the one we use, since the fully retarded amplitudes 
of the ra and R/A formalisms are the same^^. 

The four distinct types of cuts which contribute to the self-energy diagrams with two HTL 
effective vertices are depicted in Fig. ^. As just mentioned, all cut propagators are Wightman 
propagators, and they attach to the neighboring vertices like Keldysh r fields. When soft 
gluon propagators are traversed by the cut, the small P approximation G-^{P) ~ Grr{P) 
can be used. When HTL vertices are traversed by the cut, two hard propagators are put 
on-shell; as discussed in more detail in [27], the corresponding amplitudes are precisely given 
by the real-time HTL amplitudes having two external Keldysh a indices. Physically, these 
amplitudes are obtained by making the eikonal approximation in all propagators and vertices 
entering in the hard loop; the interested reader may readily verify that this reproduces the 
HTL amplitudes with two Keldysh a indices as given by the rules of section [4.1] . 

The cuts of type (i)-(iii) share a feature which is very pleasing from the viewpoint of 

^"Although the rules given by these authors deal with the imaginary part of the retarded self-energy, rather 
than the Wightman self-energy, the rule we use follows from the latter by a straightforward application of 
fluctuation-dissipation relations (KMS conditions). 

Since the most complicated self-energy diagram we need to evaluate contains "only" three loops, it is 
also possible to give a direct proof of the cutting rule in our case, starting from the standard rules of the 
Schwinger-Keldysh ra formalism applied to the calculation of the aa self-energy (which is the average of the 
two Wightman self-energies.) We have checked this; although somewhat long the proof is a succession of simple 
manipulations, which only involve the addition or subtraction of suitable closed loops of retarded propagators 
to the diagrams (such closed loops in a diagram evaluate to zero.) 
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their numerical evaluation: they are given by expressions supported on the spectral weights 
of the soft gluon propagators. In other words they split into pole-pole, pole-cut and cut-cut 
parts, according to whether the momenta Q and R are restricted to lie on the position of a 
plasmon pole, or to lie within the space-like region ("Landau cut"). This should be obvious 
from Fig. |^. In contrast, cut (iv), which represents a (two-loop!) virtual correction to the tree 
processes considered in section p^ , might be expected to induce additional complications since 
it leaves essentially unconstrained the gluon momenta that appear in it. However, somewhat 
to our surprise, under the special circumstance = we were able to bring this contribution 
into a form manifestly supported within the cut-cut region. The relevant manipulations are 
described below in greater detail. As far as we know, this additional difficulty did not show 
up in previous HTL calculations, such as Braaten and Pisarski's pioneering evaluation of 
the gluon damping rate [35], or Braaten, Pisarski and Yuan's calculation of soft dilepton 
production [36] , cut (iv) being kinematically forbidden in these cases due to the the external 
momentum being time-like. 




Figure 8: The distinct cuts which can go through the self-energy diagram (A) of Fig. [7| with two 
HTL effective vertices (drawn as loops), as explained in the text. Solid lines denote hard propagators 
and the two gluon propagators carry soft momenta. 



4.3.2 Evaluation of the cuts 

In Fig. P we give explicit expressions for the HTL effective vertices entering Fig. |8| (i)-(iii), in 
terms of certain functions, 

M-(g,i^). /f^— (4.6) 

K^%Q, P) ^ / '^^7^5{v . Q)^ , (4.7) 
J Air V ■ F 

y^'^(Q, R) = ^ [2^%/^^ + {R + PYrj''^ - (Q + PYr|^'^] . (4.8) 



which are related to the fully retarded HTL three-point vertex, its discontinuities, and to the 
tree vertex, respectively. The Keldysh indices that appear on these effective vertices, for the 
cuts (i)-(iii), are completely determined by the cutting rule we use. 

With these basic building blocks in hands, the evaluation of the cuts of type (i)-(iii) is 
relatively straightforward (our normalization is IP'^^^P) = 21mlY^{P){T /p^) with H^ as in 
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(1) 

(2) 
(3) 



JvV-P 

= m^^gf'^" X q^M^'''''{Q,R) 



v-Q- vR- 




mhgTi'f 



3 ra'bc / V V 



vP- 



2tt5{v ■ R) 



= m^gf'^" X -2K^''^'{R,P) 
mj^gf'"' X 2Ki"'{Q,P)* 



Figure 9: The HTL effective vertices that appear in Fig. ^ (i)-(iii), with ra indices shown. The 
momenta, Lorentz and color indices are as suggested by the position of these objects in Figs. ||. 



Eq. (0)): 



N,m%g^Ty2 
rr> 00 c„^ 

N,m%g^Ty2 



Q 



Q 



M^,{Q,R)M^.,.{Q,RY p^^' {Q)p^^ (-R) , 

-2M^,iQ,R)K^,,,{Q,Pr Gf {Q)p''''' i-R) 

+2K^,iQ, P)M^,,, {Q, R)* Gf {Q)p--' {-R) 
+ {Q^R), 



TT>00 



(P) 



N,m%g^Ty2 Jq 



ak,,{r,p)k;aq^p)g^r^ {Q)G7 {r) 
+{Q ^ R) , 



(4.9) 



all of which are manifestly real. We are using the approximation Grr{Q) ~ p{Q) /q^i with p = 
(Gr — Ga) being the spectral density. The contribution from cut (i) is also manifestly positive, 
as expected from its rather obvious interpretation as the square of a one-loop amplitude 
(although at this stage it might not be obvious that the sum over Lorentz indices yields a 
sum of positive terms; this is confirmed in section 4.3.3| .) In Eq. (O) the contribution from 
diagrams involving the tree interaction vertices is not shown explicitly; it can be recovered 
by the simple substitution, for each appearance of the fully retarded HTL vertex M^^{Q, R) 
or its complex conjugate. 



q^M^,{Q,R)^q'>M^,{Q,R) + V^,{Q,R). 



(4.10) 



Evaluating the cuts of type (iv) poses some additional difficulty, as mentioned above. 
The HTL diagrams contributing to it are depicted in Fig. H. The leftmost solid line in these 
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"Torg II '/ V y /'^ ryy 



"Tons II Vc/ > innr m A, , "Tnnr+P^^ IL > lonr 







Figure 10: Zoom on the ra structure of the propagators appearing in the HTL diagrams contributing 
to cut (iv) of Fig. |. 



diagrams can be identified with the cut hard loop on the left-hand side of diagram (iv) in 
Fig. |8|, collapsed to a one-dimensional line [27]: the cut (rr) HTL propagator replaces the two 
cut hard propagators in Fig. ^ (iv) , and the retarded HTL propagator on its right replaces the 
(eikonalized) third propagator of this hard loop. The reason why there appears exactly one 
cut gluon, or HTL propagator, in the rest of the diagram is because the object which stands 
on the right-hand side of the explicit cut in Fig. ^ (iv) is a one-loop retarded amplitude (in 
the HTL theory.) The direct evaluation of these HTL diagrams yields: 

nfAMiv)(^) = +2ReN,mU^T^ [ 2K^,{P,Q) 

J Q 

' M^,AQ,P)G^i^'{Q)p'"''{R)'^ 

-m,,ar,p)g%^\Q)p''''\R)'^ ^^-^^^ 

I -2K^.,,{R,P)Gf {Q)G-/ {R) + V^,,,{Q,P)Gf {Q)p-^' {R)^ . 



X < 



Using standard methods of contour integration that make use of the analyticity of the pref- 
actor K^y{P, Q) in the upper half plane, this expression can be rewritten more compactly. 
Specifically, we introduce a small ie prescription l/{r^ — ie) in all places this ap- 

pears. This is equivalent to displacing the contour of integration slightly below the real 
axis (or, equivalently, the contour of integration slightly above the real axis), and does not 
change the final answer since the numerator vanishes at = 0, p{R) being an odd function 
of r". However, the introduction of this prescription makes it possible to decompose p{R) 
into {Gr{R) — Ga{R)), the integral of each term remaining well-defined. Decomposing also 
2K{R,P) into {M{R,P) + M{—R,P)), and dropping all terms which are analytic in the 
upper-half q^ plane, one obtains: 

nfAMiv)(^) = +2ReN,m%g'T' j 2K,,{P,Q)G';t' {Q)G^/ (R) 

J Q 

^^^ f q^M,,,,{Q,P)-r^M^,,,iR,P) ^ V,>u'iQ,P) \ ^ ^^_^2) 
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Although Eq. ( [4.12| ) is valid for arbitrary , a great simplification occurs when = 0: 
the denominator l/(r*^ — ie) is antisymmetric under Q ^ R, modulo its ie prescription. 
By enforcing symmetry of the integrand under {Q ^ R) we can thus trade the prefactor 
2K{P, Q) (whose support extends to oo, a nuisance for numerical work) into the better- 
behaved combination K{P,Q) + K{P,R) = 21!ie K{Q,R), whose support lies entirely within 
the region of spacelike Q and R, \q^\ < min(g, r). Due to the explicit factor of present in the 
numerator, the ie prescription in the denominator of the HTL terms (the one involving M^/,^/) 
actually is unimportant and can be discarded, allowing {Q R) symmetry to be enforced at 
no cost. However, no such factor of q^ multiplies the term involving the tree vertex V^v'; ^-nd 
enforcing the (Q ^ R) symmetry in it gives rise to an additional contribution proportional 
to S{q'^), coming from the mismatch of ie prescriptions. Thus at p'' = Eq. ( f4.12 ) becomes: 



"y)°(iv)(^) 



8Re / [ReK^,iQ,R)]G''/{Q)G''/iR) 



M^,^,{Q,R) + 



4 V.o{q,r)GMG'^ir){S^' 



q'q^ 



lm[Kjo{p,q)-Kjo{p,r)] . (4.13) 



We remark that the apparent singularity at q^, in the term involving the tree vertex on the 
first line, is illusory (this is why we dropped the ie prescription in it.) Indeed, since the tree 
vertex with three fields is identically zero, and the tree vertex between one and two 
transverse gauge fields is explicitly proportional to q^, a singularity could only happen when 
one gluon is transverse and the other one is longitudinal. However, in that case the prefactor 
KeKiQ{Q,R) turns out to be explicitly proportional to q^. In terms of the formulae for the 
Lorentz algebra given in the next section, the corresponding statement is that Xt-l Q^- 
Thus the first line of Eq. ( [4. 131) is not sensitive to the g'^ — > region. However, although the 
prefactor Ke Ki(j{Q, R) = ^{KiQ{P,Q) + Kio{P,R)) vanishes at q^ = 0, Ki(){P,Q) itself does 
not. This is the reason why one gets a nonzero residue at g*^ = in the transverse-longitudinal 
case (and only in that case), as given on the second line of Eq. ( 4.13| ). 

To bring these expressions into a form suitable for numerical evaluation, we find it con- 
venient to decompose all individual factors into their real and imaginary parts. Doing so, the 
contributions Eq. (]49|) and the first fine of Eq. ( ^.13|) add up to (p = IReGn): 



' p{Q)p{-R) [{ReM -ReKf - (ImKQ)'^ 
-{lmKn? + {ReKf] 
+4 Im G{Q)p{R) [(Re M -ReK) Im Kq + RqK Im Kr] 
+4p(Q) Im G{R) [(Re M -ReK) Im Kr + RqK Im Kq] 
+8 Im G{Q) Im G{R) [(Re M -ReK)ReK 

-IxuKqIuiKr]. 



(4.14) 



Here we have not explicitly written the contributions involving tree vertices, which are recov- 
ered by the simple substitution Eq. ( 4.10| ), and we have dropped all Lorentz indices, which 
play no crucial role here. We are using the abbreviations M = M{Q, R), Kq = K{Q, P) and 
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Kpi = K(R,P). This expression incorporates a wealth of real and virtual physical processes, 
as discussed in Subsection |2.2| . 

In addition to Eq. ( 4.14| ) we have the contribution from the = residue, given by the 
second line of Eq. ( [4.13| ). Using Eq. ( [4.8D for Vio and performing the v integration we can 
make the latter more explicit: 



TT>00 

^^(A),gO=0 



p-q 



rrij^iP 



p^ + 



(4.15) 



Here we performed the q integration by first doing the integration over the angle between q 
and p, and then evaluating the integration over the magnitude q from its discontinuities at 
its branch cuts aX q = =bp + i[mj3,cxD). This zero-frequency contribution to Eq. (4.13) may 



appear odd-looking, compared to Eq. ( [4.251 ). However, what we regard as truly remarkable, 
is the fact that the contribution from cut (iv) (a two- loop virtual correction!), could, when 
p^ = 0, be cast into a computer-friendly form supported on the spectral weights of the gluon 



propagators. The leftover piece Eq. ( 4.15 ) seems to be the price to pay for this welcomed 
simplification. It is not clear to the authors whether such a structure persists for general 
spacelike P with p^ ^ 0. 

4.3.3 The Lorentz structure 

We have to sum over the Lorentz indices in expressions of the form 

M^,{Q,R)G^'^'' {Q)G^^\r)M;,AQ.R) ■ (4.16) 

This is where our choice of strict Coulomb gauge becomes particularly convenient. First, in 
this gauge the retarded and cut propagators have the same Lorentz structure, see Eq. (|3.4D . 
Second, in this gauge the propagator decomposes into a longitudinal part G^^ and two spatial, 
strictly transverse components, Gij = G^ (eiej+e'^e'j). We can choose one of these components, 
say, e, to lie in the plane defined by q and r and the other, e', to be orthogonal to this plane. 

These three components of the gauge propagator will give rise to four structures in 
evaluating Eq. (4.16); one contribution proportional to (G^^)"^, one contribution proportional 



to G^^G^ , and two contributions proportional to (G"^)^, one of which arises from the out of 
plane and one from the in-plane polarization states. 

The doubly longitudinal contribution to Eq. ( p^ ) is trivial; it is |M°°pG°°(Q)G°°(i?). 
Consider next the G^^ {R)G'^ {Q) contribution. To evaluate the contribution we need to study 
MiQ(Q,R), which, viewed as a vector, must involve a linear combination of and r*. Since 
the coefficient of q^ is annihilated by the Q transverse projector, only the coefficient of r* 
contributes to Eq. ( [4.16 ). We can find this contribution by applying a projector which removes 
the piece proportional to q^: 



M'^\Q, R) = r'MT-iiQ, R) + Terms proportional to q' , 



MT-LiQ,R) 



1 



2 2"(9 



rq')M' 



iO 



(4.17) 
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in which q± denotes the component of q perpendicular to p, q^p^ ^ = I'ZXrp = |rxpp. 

Now consider the contributions where both propagators are transverse. The function Mij 
vanishes when contracted against one in-plane and one out-of-plane polarization vector, by 
parity invariance in the out-of-plane direction. Therefore there are two contributions, one 
arising from the in-plane projection of Mij and one from the out-of-plane projection of Mij. 
The out-of-plane projection is 

Mt-t a{Q,R)=( - + 7, ^'^^^^ " • ^^'^^ 1 M-^- , (4.18) 

V p (11 J 

and the in-plane projection, obtained by dotting Mij against the in-plane polarization oper- 
ator for each propagator, is 

Using this procedure we find 

Eq. (|1|) = \M^^\'^G^^{Q)G^yR) + 2r\l\MT-L{Q,R)\'^G'^{Q)G^yR) 

+\MT-T,AiQ^R)\^G'^iQ)G^iR) + \Mt-tMQ,R)\^G''{Q)G^{R) . (4.20) 

We have evaluated the scalar functions entering Eq. ( [4.20| ) in terms of linear combina- 
tions of M°°(Q,i?), 1, and two new functions L{Q) and L{R), with momentum-dependent 
coefficients (that are real and analytic functions oi q^.) In a condensed notation the result 
can be written: 



Eq. (lOoD = ^PiG'^^{Q)G^^{R)Mi{Q,R)Mi{Q,Ry , (4.21) 



MiiQ, R) = X,mOO(Q, R) + YqMQ) + YR^LiR) + Zi , (4.22) 

with 



and M^^{Q,R) as defined in Eq. (|4.6| ). The sum over i in Eq. ( [4.21 ) covers the four cases 



L-L, T-L, T-T,A and T-T,B. The momentum-dependent coefficients Xi, Yq., Yh. and Z, as 
well as the prefactors Pi and choices of propagators, G*^' and G^\ are tabulated in Table 
Eq. (^). In this table we have separated the contributions to the "constant term" Z coming 



from the HTL and tree vertices. The various discontinuities of M^y which enter Eq. (4.14) 



can be obtained from the discontinuities of the basis functions entering Eq. (4.22), which are 
described in detail in Appendix 

4.3.4 Final expressions for diagram (A) 

The Wightman self-energy 11^'^'' enters the heavy quark diffusion coefficient k through: 
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Index i 


Prefactor Pj 






X, 


L-L 


1 


qOO 


(^00 


1 


T-L 








_gOp . q 


T-T,A 


1 




G^ 


1 - ql/ql 


T-T,B 


l/q\p'^q^r'^ 


G^ 


G^ 


— q^p -qp-r 



Index i 






^(HTL) 


^(troc) 


L-L 














T-L 




— q ■ r 





2/q°ml, 


T-T,A 


q^p ■ q/p^qi 


—qOp . r/p'^q'j_ 





-2/ml 


T-T,B 


1 

q q p ■ r 


n 2 
—qrp ■ q 


p'^ql 


2q • rq^p'^/ml 



Table 2: The coefncients in the expansion Eq. ( 4.21 ). The coefRcient Zi appearing in the text is the 

respectively. 



sum of its HTL and tree contributions zf^"^^^ and Z, 



(tree) 



Substituting formula Eq. ( [4.14 ) into this, upon rescaling variables by ttid and scaling out the 
prefactor from Eq. (|2.5| ) in order to obtain the dimensionless contribution to C, we obtain, 
using the decomposition Eq. ( [4.21| ) and the results of Appendix ^ 



(j4),main 



Svr 



' p^^{Q)p^^{-R) [(ReMi -Rei^O' " (Im^gJ^ 

+4 Im G«> (Q)p^> {R) [(Re Mi - Re Ki) Im Kq^ + Re Im Kr^] 
+ApQ^ (Q) Im G^'^ (R) [(Re Mi - Re Ki) Im Kr^ + Re Ki Im KqJ 
_ +8 Im (Q) Im G^^ (R) [(Re - Re /Ci) Re Ki - Im Kg, Im Kr^ 



X < 



(4.25) 



in which: 



Re Mi - Re Ki = -X 



ImA'r 



ReKi 



' p-g I \ p-r 



P\/qI-<1o 



'2q 

Y — 

^'2q 

Iu^Kr^^Yr^^^, 



-rb, — In 



q-q 



^0 



-Xi 



TT 



'^pjqi - Qo 



Yr — In f - — % ) +Zi, 



(4.26) 
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when \q^\ < q±, and: 



In 



Re Mi- Re Ki = -Xi 



2pJq^-qi 



^'2q \\q°-q\J ' 2r Vk°+^l/ 

imi.,, = x.^^^^^ + y,,f%^-,g), 

Rei^i = 0, (4.27) 



when \q^\ > q±. The arctangents on the first hne of Eq. ( 4.26D take values in [0,7r]. The 
coefficients Pi, Xi, Yq., Yr. and Zi, as well as the choices of propagator (transverse or 
longitudinal), for the different choices of i, are listed in Table |2[ Expressions for the HTL- 
resummed retarded propagators Gr and p = Gr — Ga are given in Eq. (0). The integrals 
depend only on one scale, m^,, which we have scaled out and should be set to 1 whenever 
it shows up in the formulae. As was discussed earlier, the integrand naturally splits into 
pole-pole, pole-cut and cut-cut contributions, according to whether the energies of the Q and 
R propagators lie within the space-like region (the Landau cut) or on the plasmon pole. 

The integral Eq. ( 4.25| ) is (linearly) divergent at large q, due to the transverse-transverse 



pole-pole contribution involving tree interaction vertices, which duplicates the leading-order 
gluon-scattering contribution Eq. (2.4). To obtain the correct contribution to the coefficient 



QiQCD) (defined below Eq. ( |2.5D , this leading-order contribution must be subtracted; this 
subtraction can be understood as part of a systematic matching procedure like that discussed 
in section ^ More precisely, one should subtract the contribution to Eq. ( 4.25| ) which arises 



from the tree vertices (this corresponds to keeping only the Z^^^^ part of Zi in Eq. ( 4.27 )), 
using the bare propagators. 

The \q^\ q± region presents some subtleties, that are discussed in greater detail in 
the next section: due to the square root singularities that appear in the vertex functions (in 



the terms proportional to Xi, in Eq. ( 4.2(^ )-Eq. (4.27)), which enter squared in Eq. ( [4. 25) ) 



the frequency integral is potentially logarithmically divergent in the limit q^ q±. In the 
next section we verify that the divergences cancel out between the lower and upper limits, 
\q^\ ~^ Q± and \q^\ g^, although not individually. As a consequence the integral must be 
evaluated using a Cauchy principal value prescription near q^ = q±. Actually in the next 
section we show that in addition to this Cauchy principal value integral there is an additional 



6{q — q±) type of contribution Eq. ( |4.32[ ), giving C(^A),qO=q^ — 0.023333 
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In addition to this, the zero-frequency leftover Eq. ( 4.15 ) must also be included: 



tan ^ (p) 



p2 



— ~ 0.294524 (4.28) 
32 



The evaluation of Eq. (4.25) was performed by numerical integration independently by 
the two authors. The integrals giving the cut-cut, pole-cut and pole-pole contributions are 
respectively four, three, and two-dimensional. The independent evaluations used different 
reparametrizations of the integration variables. For instance, the cut-cut integration can 
be parameterized in terms of the magnitudes of q and r, the angle between them, plus one 
frequency variable, or in terms of p, p-q, q'j_ , and one frequency variable. Both implementations 
used the Cauchy principal value prescription near = q± by "folding" the integrals in order 
that the two individually divergent parts can be added together under the integration sign and 
a convergent integral be obtained. We found satisfactory convergence in all cases, and obtain 



C(A),main ^ 0.5918. Combining with Eqs. (|42^ ) and (|3|), we thus find C(A),im ^ 0.9097 



There exist several ways to decompose Eq. ( 4.25| ) into different contributions. One way, 



although probably not gauge- fixing independent, is to separate the contributions according 
to whether they have tree-level or HTL interaction vertices. Doing so, we find that the 
contributions involving two HTL vertices are all relatively small, and add up to a relatively 
modest ~ +0.14. There are two large contributions involving two tree interaction vertices, 
both of which come from the transverse-transverse loop: one is the pole-pole contribution 
~ —0.52, which describes the influence of the plasmon dispersion relations on the scatterings, 
and the other one is the pole-cut contribution ~ +0.56, which describes scattering processes 
with the radiation or absorption of a soft plasmon. These two contributions happen to nearly 
cancel against each other, so the net contribution from diagrams with two tree interaction 
vertices is also relatively modest, ~ 0.12. The remainder of C(^) main comes from the HTL- 
tree diagrams, which add up to ~ +0.34, but originate from a large number of terms with 
different signs. A large contribution comes from the transverse-transverse cut-cut region, 
giving ~ +0.50, but this is largely cancelled by the transverse-transverse pole-cut region, 
giving ~ —0.33. Similar cancellations happen between the transverse-longitudinal cut-cut 
and pole-cut regions, which respectively give ~ —0.20 and ~ +0.18. The remainder of the 
HTL-tree contributions comes from the transverse-transverse pole-pole contribution ~ +0.30. 
We note that the total transverse-transverse pole-pole contribution, which we expect to be 
gauge invariant on its own (because this diagram is the only place where two soft transverse 
plasmons can appear) gives about ~ —0.20. 

4.3.5 A subtlelty near q^ = q± 

We now investigate in more detail the region \q^\ ~ q± of Eq. ( [4.25| ). The purpose of this 
section is to verify that the logarithmic divergences in this expression cancel out between the 
lower and upper limits, \q^\ — > and \q^\ — > q'^, so that this expression makes sense as a 
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Cauchy principal value integral. However, we will show that to take such a prescription is 



not exactly the correct thing to do, but that in addition there is the contribution Eq. (4.32). 

One procedure for regulating the divergences near \q^\ = q± is to explicitly keep the ie 
terms finite in the denominators of the HTL vertex functions M and K Eq. ( |4.6| ): in the 
time domain this regulation procedure is analogous to placing an upper bound on the time 
separation between the external legs of the self-energy diagram. At finite time separation no 
divergence is found, so that the cancellation we find in this section means that no signifi- 
cant contribution to diagram (A) arises when the time separation between the external legs 
becomes large (relative to 1/gT.) Our writing of Eq. ( [4.25| ) is entirely compatible with this 



regularization, since this expression follows from Eqs. ( f4.9D and ( 4.12 ) by simply decomposing 
each term into its real and an imaginary part; the Wightman self-energy 11^ is always purely 
real, even when this regulator is used. 

All that gets modified at finite ie are the the explicit expressions for the HTL vertex 
functions given in Eq. ( [4.26| )-( [4.27| ). The only terms we need to keep track of are those 
involving the function M^^{Q,R) and its discontinuities, e.g. the terms multiplying Xi, 
since all other terms (and propagators) are well behaved in the kinematic region \q^\ = q_\_. 
Keeping the ie's finite in the formulae of Appendix explicit expressions for the singular 
part of the various combinations of M^^{Q,R) and its discontinuities that enter Eq. ( [4.25| ) 
can be obtained: 



-P- 



Kr 



ReK 



— TT 
lip 

— IT 

4p 





(gO + ie)2 




-p • r) 




(gO - ie)2 




1 




(gO + ie)2 



+ 



+ 



+ 



9ip ■ q) 



- (gO - k)2 

9{p ■ r) 

- (gO + ie)2 

1 



(gO 



(4.29) 



Re M - Re K ~ [2e{-p ■ q) + 29{-p • r) - 1] Re . 

These expressions are valid at both positive and negative q^ , where the value they take is 
the complex conjugate of their value at positive q^ . Let us first have a look on the bracket 
from the first line of Eq. ( 4.25| ), which multiplies p{Q)p{R). Using the fact that (because 



p ■ r 



p^ —p - q) it is impossible for both p ■ q and p • r to be simultaneously negative, one can 
see from Eq. (|4.2S| ) that the singular behavior of (ReM — ReX)^ is always precisely equal to 
that of (Rei^)2. Thus the singular part of this bracket can be written: 

2(ReKQ)2 - {luiKqf - {h^Kuf = Re(/^^ + KI) (4.30) 

a result which is explicitly well-behaved and given by a Cauchy principal value integration^^ 
q_\_. The bracket on the fourth line, multiplying \m.G{Q) Im G(ii), similarly yields 

— V- + i-K5{x), which is applicable here since everything it 



near \q 



Ol 



^■^This follows from the standard formula 
multiplies is smooth as g° — > gx- 
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a finite Cauchy principal value integral. The cancellation of the divergences between the 
lower and upper limits can be seen from the explicit expressions for the divergent part of the 
bracket: 



(ReM -ReK)ReK^ ^ 2g(-p • g) + 2g(-p • r) - 1 _ 

1 ir 1 TT^ sgn (p • g)sgn • r) 
-IvuKqIthKr ~ —r 2 2 &[Q -Q±)- (4.31) 



The two lines are opposite of each other for all values of the momenta, as follows from the 
fact that p ■ q and p ■ r are never negative at the same time. 

The brackets on the second and third lines of Eq. ( 4.25| ) are manifestly finite when e = 0, 



since the real parts of M and K are only divergent for \q^\ —>■ gj, and the imaginary parts 
of the K's are only divergent for \q^\ q'^: products of these terms contain no divergence. 
However, at finite e the supports of the divergent parts of these terms overlap with each 
other, on a region of size 0{e). The contribution from this region remains finite in the e — > 
limit, giving rise to a 6{q^ — g_L)-type of contribution. It can be extracted by just taking the 
imaginary part of products of expressions from Eq. ( ^.29| ), using the formula mentioned in 
footnote upon rescaling variables by mj, one obtains the dimensionless contribution to C 
of Eq. Q: 



^' Jvi^+P) Jo P<1± , 



X < 



lp{l+p'^Y Jq p^q± 

{9{—p ■ q) — 6{—p • r)) 
X [Im G^^ {Q)p^^ {R) + p^^ (Q) Im (R)] 
+9{p ■ q)9{p ■ r) 

X [lmG««(Q)p-f^«(i?) -/>0'(Q)ImG^»(i?)] 
~ 0.0233326 , (4.32) 

which turns out to be a small contribution, most of which arising from the transverse- 
transverse contribution (this result was obtained by numerical quadrature.) 

4.4 The diagram (B) 

The diagram (B) represents the expectation value of the correlator Eq. (4.2), in which the 
heavy quark's Wilson line is expanded to linear order in the field. 



/oo rt 
dt / dt'(aM°"(t)ylO^(t')9M°^(0)), (4.33) 
oo Jo 

= ^r'' I = 0)) 27r<5(pO) , (4.34) 

3aA Jq,r r'^ + 

all fields being Keldysh r fields. The assignment of the momenta in this section is illustrated 
in Fig. 1^; P = Q + R. The fact that one of the fields carries zero frequency is probably best 
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visualized if time translation symmetry is used to move the middle field's time argument to 
t' = 0: restricting the time argument of the leftmost field to t > (thus picking up a factor 
of two), the rightmost field's time argument t" in the first line of Eq. ( 4.34 ) is then restricted 
to the range — cxd < t" < 0. However, due to the antisymmetry of the group theory factor, 
this range can be extended to cover the whole real axis, the contribution from < t" < oo 
giving zero. 




Figure 11: Zoom on the ra structure of the propagators in the HTL diagrams with topology (B), 
when the zero frequency gauge boson propagator (rightmost one) is cut (rr). These diagrams form a 
telescopic sum. 

Since there exists no bare interaction vertex with three Lorentz indices, this correlator 
only receives a contribution from the diagram with an HTL three-point vertex. The dominant 
diagrams that are allowed by the rules of the ra formalism either involve two cut (rr) gluon 
propagators and one retarded propagator, with an arr HTL vertex, or one cut propagator 
and two retarded propagators, and an aar HTL vertex. We find it convenient to organize the 
resulting diagrams into two classes, according to whether the gluon propagator which carries 
zero frequency is cut or retarded. When it is cut, one has the HTL diagrams of Fig. which 
form a telescopic sum evaluating to: 



CHmj^g'^Nc f 2ip-q 



P,Q 



1 



1 



vQ+vR- 



Cnrnj^g^N, 



Grr (p) 



2p • q 



vQ- v-R^ 
1 



(g2 + m|))(r2 + m|,) v ■ q v ■ r 



(4.35) 



The second line follows from the first by means of contour integration in the complex q^ 
plane. The fact that one ends up with an integral involving only zero-frequency propagators 
is analogous to what happened for the real part of the gluon self-energy diagrams in section 
[4.2| , which we evaluated in terms of the zero Matsubara modes. This is why we called this 
contribution "K(B)^Re"- The v integration gives the function —M^^{q, — r), which is evaluated 
in Eq. ( |A.6| ). Because it is symmetrical under q ^ r the factor 2p ■ q can be traded for a p"^. 
Upon rescaling variables by mn and scaling out the prefactor C/ff^^r^m/j/lSTr from Eq. (2.5), 
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one obtains the dimensionless contribution: 



C(S),Re = -67r 



p 



vr — cos 



1 ( 021 

qr 



P,q (1 + P'm + r'){l + q^) - (q • r)2 

0.0482933 , 



(4.36) 



in which the branch of the inverse cosine on the first hne ranges from to vr. The final result 
was obtained by means of numerical quadrature. 




4^ 




4^ 




4^ 



Figure 12: Zoom on the ra structure of the propagators in the HTL diagrams with topology (B), 
when the zero frequency gauge boson propagator (rightmost one) is retarded. Some crossed diagrams 
are not shown. 



When the zero-frequency gluon propagator in Eq. (|434| ) is a retarded propagator, one 



has the HTL diagrams of Fig. 12, the evaluation of which yields: 

'GO°(Q)G°0(i2)gOMOO(Q,i?) 
-TG?0(Q)GOO(i?)K(Q,P) 



Cnrnj^g^N, 



P,Q 



iq ■ p 



+ le q^ + te 



ir ■ p 



X < 



+TGOO(Q)GOO(i?)K(i?,P)* 
y+2TG'^^{Q)G^%R)K{P,Q). 
where we have symmetrized the factors from Eq. ( 4.34| ) under (Q R). This expression bears 



* (4.37) 



much similarity to those encountered in section 4.3 for the imaginary part of the self-energy 



diagram (A). From now on the discussion closely parallels that given there, so we will be 
brief. Indeed, the first three lines of the brace in Eq. ( [4.37 ) are very similar to the expressions 
in Eq. (|4.9| ) pertaining to cuts (i) and (ii), except that now there is only one HTL vertex 
which gets cut. These expressions are also nearly in a form suitable for numerical integration, 
since they are manifestly supported on the spectral weights of the gluon propagators (they 
decompose into pole-pole, pole-cut and cut-cut contributions, like Eq. ( [4.25|) .) However, just 



like Eq. (4.11) from section 4.3, the fourth line of Eq. (4.37) poses additional difficulty: its 



support extends beyond these regions. In section 4.3 this difficulty was dealt with by writing 
Grr{R) = {Gr{R) — GA{R))/{r^ — ie), then dropping the Ga{R) term using its complex 
analyticity in the upper-half plane, and finally using {Q R) symmetry to convert 
2K{P,Q) into the better-behaved combination {K{P,Q) +K{P,R)) = 2Re K{Q,R). When 
one repeats these manipulations here, one first needs to switch the ie prescription in the 
iq • p/{r^ + ie) prefactor in Eq. (4.37), in order to make it coherent with the one we wish to 
introduce. This gives rise to a contribution proportional to 6{q^); otherwise the results are 
very similar. 



- 29 - 



We note that Eq. ( 4.37 ) is manifestly real, the contribution from < being the complex 
conjugate of the contribution from > 0. Taking the real and imaginary parts of each term, 



as we did in section 4.3.2 , the contribution to dimensionless C Eq. ( p. 51) can be written: 

' p{Q)pi-R) [Re M{Q, R) - Re KiQ, P)] 



(B),Im 



67r 



dq^ 1 



X < 



+ 



vr 



+2p{Q) Im G{R) Im K{R, P) 
+2p(i?)ImG(Q)Imi^(Q,P) (433) 
_ +4 [Im G{Q) Im G{R) - |gO=o] Re K{Q, P) 



2pqj_ 



[Grr (q) Im G(r) + Im G(g)G„ (r)] /T . 



-0.07338 



(4.39) 



The fourth line in this expression arises solely from the fourth line of Eq. ( 4.37 ) and the sub- 
traction "Iqo" in it means to subtract Im G(Q) Im G(i?) evaluated at q^ = 0; this subtraction 
is convenient because it makes the integrand well-behaved^^ near q^ = 0, and it is justified 
by the fact that dq^K{P, Q)/{q^ + ief = 0. A part of the "leftover" from g° = on the 
last line of Eq. (|438D arises from the manipulations described in the preceding paragraph, 
concerning the fourth line of Eq. ( [4.37 ), and another part of it arises when the real part of 
the square bracket in Eq. (4.37) is taken, since this contains a 5{q^) term. 

Explicit expressions for the functions ReM, YieK and luiK that appear in Eq. ( [4.38] ) 
can be obtained by setting X = 1 and Ig = Yr = Z = in Eq. ( [42^ ) and Eq. (|]2^) 
of section 4.3.4 We evaluated Eq. ( [4.38|) by numerical quadratures; the integration splits 
into pole-pole, pole-cut and cut-cut contributions (plus the zero-frequency leftover) and their 
evaluation is very similar to the self-energy diagram (A). There is one important subtlety, 
though: neither the pole-cut, cut-cut nor q^ = Q contributions are individually well-defined, 
as they are all (logarithmically) ultraviolet divergent in the limit of large r, q fixed (taking 
r > q, for definiteness) . However their sum is UV convergent; here we skip the explicit 
verification of this fact. What this implies is that these contributions must be added under 
the integration sign, e.g. the content of the large brace in Eq. ( [4.38 ) must be added up before 
the integration over spatial q and p is performed. The largest part of Eq. ( 4.39| ) originates 
from the pole-pole region: thus diagram (B) mostly describes scatterings against longitudinal 
plasmons (this diagram is an interference term between the Coulomb and Compton channels 
for longitudinal plasmon scattering.) 

4.5 Self-energy (C): imaginary part 

The three different diagrams contributing to the cut self-energy diagram of type (C) are shown 



in Fig. 13. A direct evaluation of these diagrams would yield: 



n>00(P) = g'N^mlT / G^f,^(Q) 



2it5{v ■ P) 



27t6{v ■ {P + Q)) 



+ 



V P-v {P + QY 
2tt6{v • P) 



vP-v-P+ vP+v{P + Q)- 



(4.40) 



"'^Without this subtraction, the ie prescription in l/{q^ + ie)^ would have had to be explicitly kept. 
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Figure 13: Zoom on the ra structure of the propagators in the HTL diagrams with topology (C). 

However, this contains manifestly ill-defined factors such as 6{v ■ P)/v ■ . These ill-defined 
expressions are a typical example of pinch singularities and are due to cuts (i) and (iii) of 
Fig. 13, which are attempting to provide self-energy corrections to an on-shell propagator. 
However, as is well-known, such pinch singularities always cancel out, and in the end one 
typically obtains expression possessing a "gain-term, loss-term" type of structure character- 
istic of Boltzmann-like transport equations (see e.g. [37]). In our case the simplest way to 
regulate the pinch singularities is to use the space-time description of the four-point HTL 
with two external Keldysh a indices^^, which is simply given by an adjoint Wilson line along 
the trajectory of the light particle [27]: 

r poo rt rt' r 

U>^^{P) = -g'^N^mlT Vf,vu dte-**^-P / dt' dt" / e-^(*'-*")^-« G^;(Q). (4.41) 

Jv J-oo Jo Jo Jq 

Here t is the time separation between the two endpoints of the Wilson line and the integrand 
gives the expectation value of a Wilson line in a fluctuating background gauge field, the 
Wilson line being expanded to second order in the background field and the fluctuations 
being described by the Grr correlator. Performing the integrations over the time arguments 
and using the symmetry of Grr{Q) to drop terms which are odd under Q —Q, one obtains 
the well-defined integral: 

n> °°(P) = g'N^mlT ^ v,v^ G^^AQ) 2^ ''^" ' j^^^^T • (4-42) 

which has a transparent structure as the sum of a gain term and a loss term. 

To evaluate Eq. ( [4.42 ) we first rotate q so that its z axis lies aligned with v. The 



remaining angular integration is over the dot product n = v • p/p. Although u integration 
naturally ranges between -1 and 1, by noting that one would obtain zero if it were extended to 
cover the whole real axis, one can trade it for an integration over \u\ > 1, which we find more 
convenient for numerical purposes. Upon performing the p integration Eq. ( [4.24| ), rescaling 
variables by rrir, and factoring out a numerical prefactor, one obtains the contribution to 
dimensionless C Eq. (|2.5[) : 



^(C) = 6vr / / — / [Gf,{up + q,, (?) + (1 - ql/q^)GUw + q,,q)] /T 

Jp + P > Jl U Jg 

~ -0.132916 (4.43) 



^''The result we obtain with this regulator agrees with what one would obtain by regulating the pinch 
singularities by resumming self-energies into the propagators. Indeed, the two regularization methods only 
differ on frequency scales of order the collisional width (damping rate) F ~ g^T, whilst the effects we are 
looking at take place on the frequency scale gT. 
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where the correlators Grr are to be evaluated in the soft approximation Grr{P) = {Gr{P)- 



GA{P))T/q^. Although Eq. (4.43) can be simplified somewhat (for instance either the u or 
the p integration, at fixed up, can be done by hand), we had to evaluate it using numeri- 
cal quadrature. One encounters pole and cut types of contributions in the transverse and 
longitudinal channels. The result turns out to be almost completely determined by the con- 
tribution from the transverse pole, which produces —0.113353 by itself. Thus the important 
physics described by this diagram is not that of overlapping scattering events, but rather 
that of tree-level Coulomb scattering processes accompanied with the emission or absorption 
of a soft transverse gluon (however we do not expect diagram (C) to give a gauge-invariant 
account of these processes by itself, as the emitted transverse gluon may also be emitted at 
the exchange gluon, yielding diagrams that are included in (A).) 

4.6 The diagram (D) 

This diagram involves two insertions of a fluctuating gauge field along the heavy quark's 
adjoint Wilson line Eq. ( [4. 21 ), and can be written as: 

Hd) = o 



/ dt dt' dt"e-'*^'°e-^(*'-*")''VG'°°(P)G°°(Q). (4.44) 

Jp,Q J-oo Jo Jo 



The structure of this expression is very similar to that of Eq. ( 4.41| ), met in considering 



diagram (C), the only difference now being that the Wilson lines lie along the static trajectory 
of the heavy quark, as opposed to the light-like trajectories of the light plasma particles. Upon 
performing the time integrations, rescaling variables by itid and scaling out the numerical 
prefactor C//(7^T^m£)/187r, we obtain the dimensionless contribution: 

= evr / / ^p2 G'rUQ',P)-G^.'AO,p) ^oo^^o^ ^) 

Jp,q J ^TT 

~ 0.0675263 . (4.45) 

This was obtained by means of numerical quadrature; the integral splits into pole-pole (~ 
0.0474), pole-cut (~ 0.0261) and cut-cut (~ 0.0097) contributions, which respectively require 
one-, two-, and three-dimensional integrals. In addition there is a contribution coming from 
the zero-frequency term G^'^{0,p) (~ —0.0156), whose support in the > p region needs to 
be handled separately; these integrals posed no particular difficulty. 

This exhausts the diagrams contributing the the momentum diffusion coefficient at next- 
to-leading order in QCD. Taking the sum of the numbers we obtained we find C ~ 1.4946, 
with C as in Eq. (|2.5D . 



5. J\f=A super Yang-Mills 



Our setup for AA=4 SYM was described in section 2.4. To perform the next-to-leading order 



calculation we need the appropriate generalization of the force-force correlator and Wilson 
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line appearing in Eq. ([4.2|): 



\^ 1"^ V f 1 ab 

^{SYM)^^ (it(a'(AO + </.)"(t) Peio"^*'!^ 9'(^° + (/.)''(0)), (5.1) 

OflA J-oo L -I 



where <j) stands for the (canonically normahzed) real scalar field of SYM which couples to the 
nonrelativistic heavy quark. The derivation of this formula is identical to that of Eq. (|4.2D ; 
one must account for the infiuence of the scalar field on the eikonal propagation of the heavy 
quark and on the force it feels^^. In expanding the correlator Eq. ( |5.1| ) into powers of A9 and 
(j) only the terms with even powers of (j) have to be kept, since correlators of odd powers of 
(/)'s vanish by virtue of the global SU(4) R-symmetry of AA=4 SYM. The power-counting is 
the same as for QCD: we need to retain the diagrams with two soft loops. The nonvanishing 



ones, which we did not previously compute in dealing with QCD, are depicted in Fig. 14 



Interestingly, these diagrams are significantly easier to evaluate than those of the pre- 
ceding section, due to the simplicity of HTL amplitudes involving soft scalars: the HTL 
scalar self-energy reduces to a momentum-independent mass shift m| = AT^, and there ex- 
ists no HTL eff^ective vertex with external scalars [38]. These features are generic to theories 
containing scalar fields (with no cubic scalar self-interactions). Physically, the absence of an 
imaginary part (Landau cut) in the soft scalar two-point function (from which the momentum- 
independence of the real part of the self-energy follows, by a Kramers-Konig relation) can 
be attributed to the spin-suppression of the processes by which soft virtual scalars would 
be produced by the small-angle scattering of light plasma particles with spin. Indeed these 
processes, namely Yukawa scattering of fermions and gluon-scalar conversion, both involve a 
change in the helicity of the light particle. Since the soft virtual scalar carries no polarization 
tensor, the matrix elements for these processes must be proportional to the defiection angle, 
which is 0{g) in the HTL limit. Exactly the same mechanism prevents a background scalar 
field from interfering with the eikonal propagation of relativistic particles with spin, thus 
suppressing the insertion of external scalar legs onto existing HTL amplitudes. 

Perhaps among the most conceptually transparent sources of NLO corrections from soft 
scalars is that stemming from the imaginary part of the gluon self-energy diagram (A) of 



Fig. 14. This diagram produces a negative correction, due to the reduced phase space avail- 
able for scattering against a massive scalar. Actually, this diagram naturally combines with 
diagram (D') and with the part of diagram (B) in which the zero-frequency gluon propagator 
(which is the rightmost propagator) is a retarded propagator: their sum forms a single logical 
unit, which accounts for the effect of the scalar mass to the tree processes described by the 



^^Had we written the SYM version of Eq. ( ^.l| ) instead of that of Eq. (4.2), the force term would involve the 
manifestly gauge-covariant combination (F'" + D^4>). The spatial gauge fields A* enter this expression in just 



such a way as to give total time derivatives which can be dropped, as discussed above Eq. (4.2 
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(A) 



/ ^ / \ 

/ ^ - _ - \ 
/ \ 

/ \ 



(A') 




(B) 



(B') 





(B") 



(C) 



(D) 



(D') 



Figure 14: Additional diagrams that contribute to the next-to-leading order momentum diffusion 
coefficient in JV—i SYM. Not shown, one permutation of (D). 



second line of Eq. (|2.6D : 



A2 



247r3 Jq 



Jq/2 dEk/dk 



(g2 -I- jn2)2 



+ 



1 



1 



+ 2EI 



(5.2) 



Ek=k 



Here it is most convenient to trade the k integration for an integration over Ef^.: the integrand 
then becomes independent on the functional form of Ej^, and the massive and massless results 
only differ due to the different integration bounds at small E^. The NLO correction arises from 
the region q ~ gT, in which case the E^ integration ranging from {q/2) to y'{q/2)^ + nig can 
be done within the approximation nB{Ek){l + nB{Ek)) = T"^ / E"^. Upon rescaling variables 
by itid and factoring out A^T^mo/SGTr, one obtains the following dimensionless contribution: 




TT^ / Tdq 



-31^/2 



vV + 2-g ^ vV + 2-g 



+ 



{q" + q,/^^{q^ + 1) 

ISvr cosh-i(^/2)\ 

— ~ -0.37429 



+ 77 I 4 



1 



3 (g2 + 2)3/2 



(5.3) 



The real part of the zero- frequency self-energy diagram (A) gives zero at NLO, as can 



be seen using the argument employed in section 4.2: because the interaction vertices are 



proportional to the loop frequency, in the imaginary time formalism the diagram is saturated 
by the non-zero Matsubara frequencies, for which no 0{g) corrections arise. The tadpole 
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diagram (C) produces a (negative) momentum-independent shift to 



(5mn = 6AT 



(27r)3 



1 



1 



+ nig 



3V2 
4tt 



XTrun 



which induces a positive shift to C, 



3^/2 



(2vr)3^ -'r-rv^^27r(l+p2) 



9^/2 



0.50643 



(5.4) 



(5.5) 



A more or less related virtual contribution arises from the terms in (B) in which the 
zero- frequency gluon propagator is cut: 

A2 



6 



-2i 



P 



Grrip) Gil\Q)Gf{R) + Gf{Q)Gf,\R) 



(5.6) 



The pattern of propagators in the bracket is the same as usually arises in the calculation of 
one-loop retarded self-energies, and because the corresponding external momentum P carries 
zero frequency, one finds that the dominant result at soft p can be expressed in terms of the 
contribution from the zero Matsubara frequency. This can be derived from Eq. (^]B) by writing 



Grr{Q) = {GR{Q)—GA{Q))T/{q^—ie) in the first term, and using analyticity in the upper-half 
plane to drop the Ga{Q) contribution, thus turning this term into Gr{Q)Gr{R)T / {q^ — ie). 
Similar manipulations on the second term of the bracket yield Gr{Q)Gr{R)T/ [r^ — ie), which 
cancels against the first leaving only a (5-function at q^ = 0. Rescaling variables by mj, one 
thus obtains the dimensionless contribution: 

_ ^ r"" p^ dp 
^2^ Jo (1+F 



P 



1 



p,q 



(p2 + 1)2 (^2 + 1)(^2 + 



,2^2 



tan^ 



-1 JL 
V2 



0.37044 . 



(5.7) 

Corrections to the real part of the scalar self-energy are irrelevant at NLO, due to the 
vanishing of the imaginary part of the HTL scalar self-energy. The corrections (A') to the 
imaginary part of the scalar self-energy naturally combine with diagrams (B'), (B") and (D) to 
form a single logical unit (in diagrams (B') and (B") the zero-frequency scalar propagator must 
be a retarded propagator, again because the imaginary part of the HTL scalar self-energy 
vanishes.) Together they describe the effects of the scalar mass and plasmon dispersion 
relation on tree-level gluon-scalar scattering, as well as new processes, involving physical 
longitudinal plasmons in the external states or overlapping scattering events (associated with 
gluon propagators in the Landau cut.) Processes involving longitudinal gluons also occur 
in a Compton-like channel, which interferes with the scalar-exchange channel: this is what 
diagrams (B'), (B") and (D) account for. Using formulae Eq. ( 4.34| ) and Eq. ( [4.45| ) these 
diagrams add up to: 



p,L 



~ -0.31086. 



^p'GUQ) 



P 



(p-qy 



+ p'gOO(Q) 



(p2 + 1)2 



p2 + 



^1 -l.GOO(g,0) 



(5.8) 
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The contribution from the transverse gauge field is hnearly divergent at large momenta, 
where it degenerates to the tree- level contribution which was already included in Eq. (|2.6D ; 
this must be subtracted. What we have integrated numerically is the difference between 
Eq. ( |5.8| ) and the same expression evaluated with the tree Grr propagators, which is given by 
a convergent integral. The integral splits into pole-pole and pole-cut contributions, with the 
largest contribution arising from the transverse gluon pole, which gives ~ —0.283. 

This exhausts the list of NLO contributions whose origin is proper to the SYM theory. 
Taking the sum of Eqs. 5.3, 5.5, 5.7 and 5.5 we obtain the next-to-leading order coefficient 
fjisYM) ^ 0.19172. 
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A. The functions M°°{Q,R) and L{Q) 

The functions M^^(Q, R) and L{Q) have appeared in some form or the other in previous work 
( [35] [36] [18], to cite a few.) The function L{Q) is familiar from the longitudinal HTL gluon 
self-energy: 



HQ) 



dQy 1 



Air V ■ Q 

- m 

2q \q^' — q + ie 



^iM^o^^). (A.l) 



The branch of the logarithm is such that this function is real for time- like Q, and acquires a 
positive imaginary part at space-like Q, the so-called Landau cut. 

To understand the analytic structure of M^^{Q,R) we find convenient to express it as a 
sum of two terms, 

M^'^{Q,R) = M^^{P,Q) + M^^{P,R), (A.2) 



and to analyze those two terms separately. Actually, the splitting of Eq. (A.2) as a sum 



of two terms is precisely how the function M'^'^{Q, R) arose in the first place, in the HTL 
effective vertices of Fig. ^, so in a sense it is rather natural. In [18] an expression for the 
(Lorentz-invariant) function M^^(P,Q) was obtained by first combining the denominators 
1/v ■ P~ and 1/v ■ by means of standard Feynman parameterization: 

^ f 1 -1 
dx / TT = / dx- 



Jv [xv ■ Q + {I — x)v • P — ie]^ Jo {xQ + (1 — x)P)^ — ieq^ 

(A.4) 



-1 ^^( Q-P + V^ 
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where: 



Q2p2 _ . py 



(A.5) 



Proper care must be given to the ie prescriptions in Eq. (A.4), since we are interested in 
this function and its discontinuities at physical values of the momenta. The correct procedure 
follows from noting that, from its definition, the function M^^{Q,P) is analytic in the upper 
half complex plane. Its integral representation is unambiguous for Q time-like, in which 
case it has no discontinuity across the real function of q'^, with fixed p and 

q. Using a change of variables v — v in the integral representation Eq. ( A. 3 ), one shows 
that flipping the sign of in M^^{Q,P) is equivalent to complex conjugation. Explicit 
expressions, when g'^ > 0, are: 



-1 



2pV9o -91 
-1 



-1 



In 



In 



p-q + P\Jql-q\ 



-p-q+P\lql -q\^ 



tail 



p-q + P\Jql - q\ 



^p-q-p\/qo -ql 



.1 / P\/<il-<lo 

p-q 



+ in 



+2TTi9{-p ■ q) 



> 



(A.6) 



< < q±, 



where all logarithms are real, and the arctangents range from to vr. 

The equations Eq. ( |A.6| ) can be understood from Eq. ( [A.4| ) as follows. When q^ > q, the 



magnitude of the square root is automatically larger than |g -p], thus the denominator in the 



logarithm in Eq. (A.4) is negative. Since this denominator has a small negative imaginary 
part, the logarithm acquires a positive phase +Z7r. This choice of branch can be verified from 
the large q^ limit, where M^^{Q, R) — > —iixj^pcP according to its definition Eq. ( |A.3| ). When 
cP crosses q {Q becomes spacelike), the square root becomes equal to |g • p|, hence either 
the numerator or the denominator of the logarithm vanishes. Which one vanishes depends 
on the sign of q ■ p, explaining the appearance of the 0{—q ■ p) function in the second case. 
Finally, when q^ < q±, the square root becomes imaginary and the logarithm goes over to an 
arctangent. 

The function K^^{Q,P) is (half) the discontinuity of M°°(Q,P) across the real o 
axis, and can be extracted from the previous results by writing K{Q, P) = {M{Q, P) + 
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M(-g,P))/2: 




q<q°, 

9± < g° < q, 



(A.7) 



— TT 



< g° < q±. 
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